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PREFACE 

Many  numerical  models  used  In  meteorology  and  other  branches  of 
geophysics  produce  large  numbers  of  fields  of  data  that  can  best  be 
studied  by  constructing  graphical  contours  of  the  fields.  This  Is 
particularly  true  of  the  cumulus  dynamics  program  developed  at  RAND 
under  sponsorship  of  the  Naval  Research  Laboratory.  A  given  computer 
run  of  this  program  may  produce  fields  of  five  or  more  variables  at 
upwards  of  100  time  steps.  If  any  substantial  fraction  of  the  fields 
Is  to  be  contoured,  It  must  be  done  mechanically. 

Contour  programs  have  been  written  for  use  In  most  numerical  weather 
prediction  facilities  as  well  as  In  some  research  agencies.  Unfortunate¬ 
ly,  they  are  all  specialized  as  regards  either  the  fields  contoured  or 
the  equipment  on  which  the  contouring  was  done,  or  both.  When  attempts 
to  find  a  ready-made  contour  program  suitable  for  the  RAND  cumulus 
dynamics  model  failed,  It  was  decided  to  develop  one  along  the  simplest 
possible  lines  that  would  give  reasonably  accurate  and  esthetically 
pleasing  results.  The  method  so  developed  Is  described  herein. 
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ABSTRACT 


Problems  of  graphical  display  of  a  dependent  variable  aa  a  function 
of  two  Independent  variables  are  discussed.  A  procedure  for  objective 
construction  of  contours  is  described,  and  the  program  logic  Is  presented. 
Examples  of  contours  of  fields  developed  by  a  cumulus  dynamics  program 
and  produced  by  a  General  Dynamics  S-C4060  are  shown. 
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I .  INTRODUCTION 

In  many  geophysical  problems  Che  values  of  a  dependent  variable 
are  best  portrayed  as  a  function  of  two  Independent  variables.  This 
is  most  conveniently  done  by  drawing  contours  of  the  dependent  variable 
on  a  grid  representing  the  Independent  variables.  Ordinarily  the  value 
of  the  dependent  variable  Is  known  for  only  a  limited  number  of  pairs, 
and  so  some  form  of  Interpolation  Is  necessary.  The  most  direct  pro¬ 
cedure  Is  to  plot  the  known  values  on  a  grid,  and  then  to  draw  the  con¬ 
tours  by  hand,  Interpolating  by  eye.  Inevitably  some  smoothing  results, 
the  amount  being  determined  subjectively  and  often  varying  from  one  part 
of  the  chart  to  another.  After  skill  has  been  acquired  by  considerable 
practice,  this  procedure  can  result  in  an  esthetically  pleasing  chart 
that  will  yield  the  value  of  the  dependent  variable  at  any  point  In 
the  field  to  a  degree  of  accuracy  consistent  with  the  accuracy  and 
density  of  the  original  data. 

This  procedure,  however,  is  time-consuming,  and  moreover  requires 
the  acquisition  of  a  specialized  skill.  When  the  data  concerned  are 
the  output  of  a  computer  program,  there  are  frequently  many  charts  to 
plot  and  analyze,  and  the  Investigator  is  often  deterred  by  the  sheer 
magnitude  of  the  task.  Some  help  may  be  afforded  by  programming  the 
computer  to  print  the  values  of  the  dependent  variable  In  their  appro¬ 
priate  spatial  relationships,  thereby  eliminating  the  tedious  step  of 
extracting  values  from  a  table  and  plotting  them  on  a  grid.  This  pro¬ 
cedure,  however,  does  nothing  to  reduce  the  more  difficult  and  specialized 
task  of  analysis.  In  some  Instances  the  computer  has  In  addition  been 
allowed  to  Interpolate  among  the  grid  values  and  the  printer  has  been 
programmed  to  strike  characters  in  the  appropriate  location,  so  that 
the  analyst  need  only  connect  like  characters  to  draw  the  contours. 

The  results,  however,  are  at  best  crude  and  inelegant. 

A  better  method  is  to  use  a  specialized  graphical-output  device. 

Two  categories  of  these  devices  are  in  general  use:  mechanical  and 
electronic.  With  the  former  the  computer  determines  by  interpolation 
the  coordinates  of  successive  points  on  a  contour  at  close  intervals, 
and  the  output  device  thereupon  causes  a  pen  to  move  accordingly  across 
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a  sheet  of  paper.  With  the  latter  the  computer  also  determines  the  co¬ 
ordinates  of  successive  points,  but  the  contour  segment  Is  then  displayed 
electronically  on  a  screen  that  may  be  photographed  for  permanent  record. 
Both  usually  have  means  of  producing  labels  and  other  alphanumeric 
Information.  The  essential  difference  from  the  point  of  view  of  the 
programmer  Is  that  with  the  mechanical  curve -follower  It  Is  necessary 
chat  all  the  points  of  a  given  contour  be  sequentially  determined  before 
the  next  contour  is  considered,  whereas  the  electronic  device  permits 
the  determination  of  all  contour  segments  within  a  given  area  (say  one 
grid  Bquare)  before  the  next  area  is  considered.  The  latter  system 
leads  to  significantly  fewer  logical  complexities.  On  the  other  hand, 
the  Inertia  of  the  pen  arm  of  a  mechanical  device  may  cause  the  contour 
to  be  smoothly  curved  (which  Is  in  accordance  with  our  notion  of  how 
many  variables  behave  in  nature),  whereas  the  electronic  device  usually 
displays  a  series  of  straight-line  segments  with  angular  joinings,  unless 
one  resorts  to  complicated  programming.  In  the  present  study  the  General 
Dynamics  S-C4060,  an  electronic  device,  was  used. 
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II .  INTERPOLATION 

If  Che  dependent  variable  Is  conceived  Co  be  a  smooch  function  of 
Che  independent  vcriables,  some  form  of  surface  fitting  is  indicated 
for  interpolation  among  grid  points.  A  quadratic  surface  can  be  fitted 
exactly  to  six  grid  points,  but  since  it  is  impossible  to  select  six 
points  symmetrically  distributed  about  a  grid  square,  it  is  better  to 
use  twelve  or  sixteen  points  (see  Fig.  1)  and  to  fit  the  surface  by 
the  method  of  least  squares.  Of  course  with  these  many  grid  points  a 
higher-order  grid  surface  could  be  fitted,  but  the  quadratic  surface  is 
likely  to  give  the  smoothest  curves,  and  without  further  knowledge  of 
the  nature  of  the  dependent  variable,  no  particular  surface  can  be 
definitely  said  to  be  better  than  another.  One  drawback  of  arrays 
such  as  those  of  Fig.  1  is  that  they  do  not  lend  themselves  to  interpola¬ 
tion  within  grid  squares  adjacent  to  the  boundaries.  Another  disadvan¬ 
tage  of  the  surface-fitting  method  is  that  although  the  resulting  curves 
are  smooth  within  the  central  grid  square,  there  is  no  assurance  that 
they  will  Join  smoothly  (or  even  join  at  all)  with  curves  in  an  adjacent 
grid  square,  which  are  computed  from  a  different  set  of  grid  points. 

If  we  are  willing  to  accept  angular  joinings  at  grid  lir.-»s  and 
straight  contour  segments  between,  we  can  effect  a  great  simplification 
by  using  linear  interpolation  along  the  edges  rather  than  surface  fitting. 
If  the  grid  mesh  is  sufficiently  fine  and  the  dependent  variable  is 
reasonably  well-behaved,  it  is  possible  to  generate  quite  acceptable 
contours  by  this  method,  even  though  the  finished  chart  is  less  pleasing 
in  appearance  than  one  drawn  by  e  skillful  analyst. 

The  method  chosen,  therefore,  is  to  scan  all  the  grid  squares  once 
for  each  value  corresponding  to  a  contour  to  be  drawn,  to  find  by  linear 
interpolation  all  the  points  where  the  contour  intersects  the  edge  of 
the  grid  square,  and  to  instruct  the  graphical-output  device  to  connect 
these  points  with  straight-line  segments.  The  procedure  for  doing  this 
and  the  resolution  of  certain  ambiguities  that  arise  are  described  below. 

*Double  linear  Interpolation  carried  out  in  the  interior  of  a  grid 
square  is  equivalent  to  fitting  a  quadratic  surface  without  the  cross- 
product  term.  The  resulting  contour  segments  are  portions  of  hyperbolas. 
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Fig.  1  —  Grid-point  arrays  suitable  £or  interpolation  within  shaded 
square.  Solid  dots  are  Cor  12-point  array;  circles  are 
added  for  16-polnt  array. 
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III.  ENUMERATION  OF  POSSIBLE  CASES 
Let  the  following  symbols  have  the  meanings  shown: 

G  ■  value  of  the  dependent  variable. 

C  ■  value  of  the  contour  under  consideration, 
i  or  j  ■  subscript  designating  a  grid  point.  By  analogy  with 
a  map,  the  subscript  has  the  value  1,  2,  3,  or  4,  re¬ 
spectively,  for  the  "northwest,"  "southwest,"  "north¬ 
east,"  and  "southeast"  corners  of  a  grid  square. 

(See  numbers  associated  with  the  shaded  square  in 
Fig.  1.) 


In  general,  if  6^  *  0,  the  contour  will  pass  through  the  grid 
point  designated  by  i.  (Exceptions  to  this  rule  will  be  discussed 
later.)  If  1  and  j  represent  two  grid  points  on  one  edge  of  a  grid 
square,  and  6^  and  6j  have  opposite  signs,  the  contour  will  Intersect 
that  edge.  Linear  interpolation  places  the  point  of  intersection  at 
a  distance  (C  -  -  G^) ,  in  units  of  mesh  length,  from  the  point 

i.  Continuity  requires  that  if  one  edge  is  so  Intersected,  the  con¬ 
tour  must  also  Intersect  one  of  the  other  three  edges  or  else  go  through 
one  of  the  two  grid  points  not  on  the  first  edge.  The  contour  is 
assumed  to  follow  a  straight  line  between  the  two  endpoints  so  determined. 
The  four  values  6^,  6j,  6^,  snd  6^  give  all  the  information  that  is 
available  or  (if  certain  conventions  are  adopted)  necessary  to  define 
the  course  of  the  contour  within  the  grid  square. 

Much  useful  information  is  given  merely  by  the  signs  of  6^,  6^, 

6^,  and  6^.  The  type  of  contour  configuration  within  a  grid  square, 
though  not  necessarily  its  exact  description,  is  given  uniquely  by 
these  four  bits  (actually  ternary  digits)  of  information.  The  number 
of  order-preserving  combinations  of  m  members  in  which  each  member  is  . 
one  of  n  symbols  is  one  greater  than  the  largest  m-digit  number  to  number 
base  n;  i.e.,  it  is  nm.  Since  in  the  present  instance  m  ■  4,  representing 
the  four  corners  of  the  grid  square,  and  n  -  3,  representing  the  symbols  -, 
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0,  and  +,  the  number  of  combinations  is  81.  Since  (C  -  G1>/(G  -  G^)  - 
6 ^ 1  "  6j)>  it  is  obvious  that  if  all  the  signs  of  the  quantities 
®1»  ®2*  ®3*  an<*  ^4  are  reversed,  the  course  of  the  contour  within  the 
grid  square  remains  unchanged.  Thus  about  half  of  the  81  cases  are 
redundant  and  need  not  be  considered  separately.  The  remaining  41  cases 
are  illustrated  in  Fig.  2. 

In  this  figure  the  sign  of  6^  is  shown  beside  the  grid  point.  The 
square  in  row  A  and  column  e  (designated  Ae) ,  having  each  6^  equal  to 
zero,  is  unaffected  by  a  change  of  sign.  If  A6  were  illustrated,  it 
would  be  the  negative  of  A£,  Ay  would  be  the  negative  of  Af),  and  so  on. 

The  40  omitted  cases  are  symmetrical  with  the  40  nonzero  cases  illustrated 

The  contours  are  shown  schematically  in  Fig.  2  as  heavy  lines.  When 
a  contour  Intersects  an  edge,  it  is  shown  in  this  Illustration  as  inter¬ 
secting  the  midpoint,  although  in  practice  its  intersection  would  be 
determined  by  the  magnitude  of  6^/(6^  "  6^).  In  four  cases  (All,  B6 , 

Cs,  and  D0)  the  ratio  6^/<6 ^  -  6j)  has  the  form  0/0,  but  it  is  arbi¬ 
trarily  assumed  to  take  the  value  1/2.  Such  ambiguities  are  di tcuased 
more  fully  in  the  next  chapter. 

In  order  to  construct  contours  as  illustrated  in  Fig.  2,  the  pro¬ 
gram  must  compute  raster  numbers  for  the  two  coordinates  of  each  end¬ 
point  of  each  contour  segment.  For  the  x-coordinates  there  are  five 
conditions,  including  the  arbitrary  case  just  mentioned,  each  requiring 
its  own  method  of  computation.  They  are: 


1. 

2. 

3. 

4. 

5. 


The  endpoint  lies  on  the  left-hand  edge. 

The  endpoint  is  determined  by  6^/(6^  ~  63)* 

The  endpoint  is  determined  by  63^2  *  ®4^’ 

The  endpoint  lies  on  the  right-hand  edge. 

The  endpoint  lies  midway  on  the  upper  or  lower  edge. 

The  corresponding  conditions  for  the  y-coordlnate  are: 

1.  The  endpoint  lies  on  the  upper  edge. 

The  endpoint  is  determined  by  6^/(6^  ~  ftj) • 

The  endpoint  is  determined  by  6 3/(63  -  6^). 

The  endpoint  lies  on  the  lower  edge. 

The  endpoint  lies  mldwsy  on  the  right-hand  or  left-hand 
edge. 


2. 

3. 

4. 

5. 


Fig.  2  Possible  contour  configurations  in  a  grid  square. 


-8- 


Although  for  both  coordinates  conditions  1  and  4  are  special  cases 
of  2  and  3,  it  is  desirable  to  consider  them  separately.  Thus  the 
general  nature  of  the  contour  configuration  can  be  concisely  expressed 
by  a  four-digit  number,  1^  Ij  where  1^  and  J^,  each  running  from 
1  to  5,  refer  respectively  to  the  x-  and  y-coordlnates  of  one  endpoint 
and  I 2  and  similarly  represent  the  other  endpoint.  This  representa¬ 
tion  corresponds  exactly  to  the  schematic  graphical  representation,  but 
is  more  convenient  to  use  in  blocking  out  the  program.  The  digital 
representation  is  shown  in  Fig.  2  below  the  grid  squares.  Some  cases 
have  two  contour  segments  and  some  have  none.  If  there  are  none,  the 
graphical  representation  is  a  zero  in  the  grid  square,  and  the  digital 
representation  is  0000. 


IV.  RESOLUTION  OP  AMBIGUITIES 


In  many  instances  (e.g.,  BP,  0»,  CTU  Dar)  the  placement  of  the  con¬ 
tour  is  unequivocal.  In  other  cases,  there  are  several  logical  possi¬ 
bilities,  and  the  choice  made  was  more  or  less  arbitrary.  Occasionally 
information  from  adjoining  grid  squares  is  needed  to  make  the  best 
choice,  but  it  vaa  felt  that  the  improvement  in  pattern  gained  by  using 
such  Information  is  not  worth  the  extra  effort  entailed. 

The  most  ambiguous  case  is  that  of  the  saddle  point,  Cy.  On  the 
basis  of  the  available  data,  i.e.,  the  value  of  G  at  the  four  corner 
points,  any  of  the  three  configurations  shown  in  Fig.  3  is  possible. 

If  the  contours  intersect  the  edges  near  their  midpoints,  there  Is  no 
logical  basis  for  choice  between  configurations  I  and  II,  leaving  III 
as  the  best  choice,  except  for  the  unaccountable  aversion  many  analysts 
have  to  drawing  crossed  contours .  (It  can  be  shown  that  for  every  saddle 
point  there  exists  a  contour  value  that  must  have  an  intersection.)  If, 
however,  the  contours  intersect  the  edges  near  one  or  more  grid  points, 
configuration  III  Is  not  necessarily  the  most  logical  choice,  and  there 
may  be  a  rational  basis  for  a  choice  between  configurations  I  and  II. 

The  various  possibilities  can  all  be  tested  by  appropriate  programming, 
and  suitable  choices  can  be  made  objectively,  leaving  but  a  small 
residual  ambiguity;  however,  In  view  of  the  complexities  introduced  by 
this  procedure  and  the  relative  Infrequency  of  occurrence  of  saddle 
points,  It  was  felt  that  such  programing  was  not  justified.  According¬ 
ly,  a  completely  arbitrary  decision  was  made:  configuration  III  was 
chosen  for  all  saddle  points. 

Other  ambiguous  cases  all  involve  a  value  of  zero  for  6^.  In  cases 
such  as  BT|  the  value  of  zero  was  ignored  on  the  assumption  that  if  it 
were  of  any  consequence  it  would  be  better  taken  into  account  in  an 
adjacent  grid  square.  Cases  such  as  B6  occur  mainly  on  the  border. 

In  this  configuration  the  endpoint  of  the  interior  contour  segment 
could  logically  lie  anywhere  on  tue  zero  edge,  but  without  knowledge 
of  the  course  of  the  contour  segment  connecting  with  the  other  endpoint 
and  lying  in  the  adjoining  grid  square,  the  most  reasonable  choice  of 
a  terminator  is  the  midpoint  of  the  zero  edge. 


Some  fields  have  extensive  flat  regions  (G  ■  constant),  and  fre¬ 
quently  C  is  selected  so  that  C  “  G  in  such  a  region.  A  common  example 
is  a  field  in  which  G  is  nonzero  in  only  a  limited  region,  in  which 

case  C  *  0  is  a  desirable  choice  for  a  base  contour.  Case  Ae  is  re¬ 

peatedly  found  in  the  midst  of  a  flat  region.  Obviously  in  this  circum¬ 
stance  a  contour  should  not  be  drawn  along  the  edges  of  each  grid  square. 
Cases  such  as  AC,  B£,  and  B6  are  found  near  the  boundaries  of  a  flat 

region.  Thus  a  flat  region  having  the  value  C  will  be  surrounded  by  a 

contour  of  the  same  value. 


Fig.  3  —  Possible  contour  configurations  for  saddle  point 


V.  PROCEDURE 


In  practice  the  entire  grid  is  scanned  once  for  each  contour-value 
desired.  For  each  grid  square  the  signs  of  6^,  6g»  6^,  and  6^  are 
determined.  If  6^  >  0,  the  appropriate  case  is  selected  from  rows  C, 

D,  and  E  of  Fig.  2.  If  6^  <  0,  the  sign  of  each  6^  is  changed,  and 
again  a  selection  is  made  from  rows  C,  D,  and  E.  If,  however,  6^  •  0, 
the  sign  of  &2  is  examined.  If  >  0,  a  choice  is  made  from  row  B; 
if  <  0,  all  signs  are  changed,  and  a  choice  is  made  from  row  B; 
but  if  both  6^  ■  0  and  62  ■  0,  6^  is  examined.  For  6^  1*  0,  one  of  the 
cases  A'H,  A0,  or  A(  is  chosen,  and  signs  are  changed  if  necessary. 
Finally,  if  6^  ■  62  *  6^  “  0,  case  Ac  is  chosen  for  6^  ■  0  and  A£  for 
6^  /  0.  Once  the  applicable  case  has  been  chosen,  it  is  a  simple 
matter  to  determine  in  terms  of  raster  number  the  coordinates  of  the 
two  endpoints  of  each  contour  segment  in  the  grid  square  and  to  require 
the  display  device  to  construct  a  straight  line  between  palrB  of  coor¬ 
dinates  . 

A  flow  diagram  of  the  procedure  is  given  in  Fig.  4.  The  initial 
elimination  of  cases  carried  out  by  adjusting  the  signs  of  the  6^  is 
shown  in  Fig.  4a,  and  the  specification  of  configuration  is  shown  in 
the  remaining  parts  of  Fig.  4.  In  this  figure  a  four-character  symbol 
in  a  rectangle  indicates  a  procedure  for  finding  raster  numbers  for 
endpoints  of  a  contour  segment.  A  digit  from  1  to  4  in  a  particular 
position  indicates  that  the  raster  number  of  the  coordinate  correspond¬ 
ing  to  that  position  is  to  be  found  in  accordance  with  the  digital  code 
of  Fig.  2.  An  X  indicates  that  no  new  raster  number  is  computed;  the 
one  previously  determined  for  that  position  stands.  The  letter  L  in 
a  rectangle  indicates  that  the  display  device  is  directed  to  draw  a  line 
between  the  points  (1^,  J^)  and  (Ij,  Jj)  currently  specified.  The  com¬ 
plete  process  described  by  Fig.  4  is  performed  for  each  contour  value 
prescribed. 

Since  several  configurations  have  contours  along  the  edges  of  the 
grid  squares,  it  is  conceivable  that  some  contour  segments  might  be 
drawn  twice.  For  example,  case  A9  could  lie  immediately  above  case  Be. 
In  this  eventuality  the  finished  contour  will  merely  have  a  segment 


Specify  C 


Fig.  4b  —  Flow  diagram. 


Fig.  4c  —  Flow  diagram. 
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Flg.  4e  --  Flow  diagram. 
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darker  than  the  rest  of  it. 

Among  the  many  possible  juxtapositions  of  the  41  configurations 
of  Fig.  2,  some  will  occasionally  produce  contours  different  from 
those  that  would  be  hand-drawn  by  an  analyst.  One  possibility  is  shown 
in  Fig.  5.  The  analyst  would  probably  draw  some  kind  of  nose  or  loop 
rather  than  the  "dangling"  contour  between  the  squares  labeled  Demand 
Be,  its  exact  shape  depending  on  the  magnitudes  of  6^  in  the  De  square 
and  &2in  the  Be  square.  For  exrmple,  if  the  latter  is  near  zero,  the 
analyst  might  draw  a  curved  contour  segment  (shown  as  a  dashed  line) 
in  place  of  the  segment  common  to  the  Be  and  Ad  squares.  In  this  case 
the  objective  program  is  concerned  only  with  whether  6^  of  the  Be 
square  is  zero  or  nonzero,  and  ignores  its  magnitude.  Numerous  other 
cases  could  be  constructed  in  which  the  objective  contours  have  un¬ 
desirable  features,  but  in  practice  their  occurrence  will  probably  be 
Infrequent. 
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Fig.  5  —  Example  of  a  "Jangling"  contour. 
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VI.  ORIENTATION  AND  LABELS 

Unless  they  are  oriented  with  respect  to  a  grid  and  Identified  In 
some  way,  the  contours  are  of  little  value.  The  grid  ordinarily  con¬ 
sists  of  full  lines  entered  at  some  Interval  prescribed  by  the  program¬ 
mer  and  labeled  in  the  margin.  If  that,  set  of  line?,  together  with  the 
contours,  causes  too  much  clutter,  only  the  Intersections  of  the  grid 
lines  need  be  designated,  or,  simplest  of  all,  mere  tick  marks  along  the 
boundaries  will  sometimes  suffice. 

Most  of  the  devices  have  means  whereby  numbers,  letters,  or  other 
symbols  can  be  displayed.  ThiB  is  Important,  as  labels  are  essential. 
Among  the  labels  required  for  almost  any  application  are  the  following: 

a.  Title  or  other  identification  of  the  whole  chart, 

Including  pertinent  information  such  as  date  and  time. 

b.  Grid  labels  sufficient  to  identify  the  coordinates  of 
any  point  on  the  chart. 

c.  Indication  of  the  value  of  G  at  significant  points, 
such  as  maxima  and  minima. 

d.  Indication  of  the  value  of  the  contours  shown. 

Few  difficulties  are  encountered  in  handling  the  first  two  items. 

The  other  two,  however,  are  less  simple.  For  example,  how  are  the 
extrema  of  G  to  be  defined?  In  the  first  place,  we  can  distinguish 
between  an  absolute  and  a  relative  extremum.  An  absolute  maximum, 
for  example,  could  be  a  value  that  is  not  exceeded  at  any  point  in 
the  domain,  whereas  a  relative  maximum  could  be  a  value  that  is  not 
exceeded  at  any  point  in  some  neighborhood  of  the  point  in  question. 

By  these  definitions  the  extrema  are  not  unique;  l.e.,  the  extreme 
value  can  be  shared  by  several- adjacent  points.  A  good  example  of 
this  is  a  field  that  is  positive  in  some  limited  region  and  zero  else¬ 
where.  In  the  large  flat  region  of  zeroes,  each  point  is  an  absolute 
minimum,  yet  it  would  be  most  undesirable  to  label  each  as  such. 

The  way  chosen  to  solve  this  problem  was  to  define  as  ■  relative 
maximum  a  grid-point  value  of  G  that  exceeds  the  value  at  each  of  the 
eight  nearest  grid  points  (with  the  converse  definition  for  a  relative 
minimum),  and  to  label  only  those  absolute  extrema  that  are  also  relative 
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extrema.  To  be  sure,  some  bonaflde  relative  extrema  will  be  overlooked 
by  this  process  (as  where  two  adjacent  grid  points  have  the  same  extreme 
value),  but  this  Is  thought  to  be  preferable  to  the  labeling  of  several 
adjacent  points  in  a  flat  region. 

Even  so.  If  G  does  not  vary  smoothly,  there  may  be  many  inconse¬ 
quential  relative  extrema  selected  in  a  nearly  flat  region.  It  may 
not  be  desirable  to  reduce  the  incidence  of  these  extrema  by  smoothing 
the  field  of  G  before  contouring,  yet  the  inclusion  of  the  location 
and  value  of  each  of  them  would  clutter  the  chart.  One  possible  com¬ 
promise  is  to  Indicate  the  location  of  each  relative  extremum  with  an 
appropriate  symbol  (as  GD  for  a  maximum  and  0  for  a  minimum) ,  but  to 
give  value  labels  only  to  those  relative  extrema  that  are  also  absolute 
extrema  or  that  differ  in  value  from  the  nearest  contour  by  more  than 
a  specified  amount. 

Finally,  the  value  of  the  contours  should  be  indicated.  Ideally, 
each  contour  should  be  labeled  at  one  point  and  at  a  few  additional 
points  if  it  is  long  or  complicated  in  shape.  A  program  that  follows 
each  contour  from  beginning  to  end  can  do  this  easily  enough,  but  one 
that  treats  each  contour  as  a  number  of  short  segments  rather  than  as 
an  entity  cannot.  It  would  be  intolerable  to  label  each  segment,  and 
complex  and  cumbersome  to  program  logic  to  determine  which  one  of  the 
many  segments  constituting  a  contour  should  be  labeled.  Fortunately, 
the  labeling  of  contours,  although  highly  desirable,  is  not  mandatory. 
If  the  values  of  a  base  contour  and  of  the  uniform  contour  interval 
are  specified  on  the  chart,  together  with  the  locations  and  values  of 
Important  extrema,  the  value  of  any  particular  contour  can  be 
Inferred  by  counting. 
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VII .  EXAMPLES 


The  procedure  described  herein  has  been  programmed  for  use  with 
the  General  Dynamics  S-C4060.  Some  examples  of  fields  computed  by  a 
cumulus  dynamics  program  are  shown  In  Figs.  6,  7,  8,  and  9. 

At  the  top  of  each  chart  la  a  line  describing  the  field  and  giving 
the  simulated  time  or  "cloud  time"  In  minutes  from  the  start  of  computa¬ 
tion.  On  the  same  line  the  contour  Interval  la  shown.  Although  a  base 
contour  value  and  the  contour  Interval  are  specified  by  the  programmer, 
the  program  may  multiply  the  specified  interval  by  a  power  of  2  in  order 
to  avoid  having  too  many  or  too  few  contours . 

The  grid  is  so  designed  that  regardless  of  the  mesh  length  of 
the  given  data,  lines  are  printed  and  labeled  at  Intervals  of  1000 
units.  Additional  boundary  lines  are  also  printed  as  needed.  For 
example,  in  the  data  as  presented  to  the  contour  program  in  the  present 
example,  grid  points  were  at  abscissae  of  -100,  100,  300,  ...,  5500, 
5700.  Printed  lines  are  at  -100,  0,  1000,  ....  5000,  5700. 

Maxima  are  indicated  by  the  symbol  [+) ;  minima  by  ©.  They  are 
labeled  only  if  they  are  absolute  extrema  or  If  they  differ  from  the 
nearest  contour  by  more  than  a  tenth  of  the  contour  Interval.  Contours 
are  not  labeled. 

Figure  6  shows  a  slightly  "rough"  field;  that  is,  one  with  numer¬ 
ous  insignificant  extrema  and  some  contour  irregularities.  In  the  upper 
right-hand  part  of  the  chart  are  a  "dangling"  contour,  some  heavy, 
doubly  struck  segments,  and  other  oddities.  They  occur  in  a  relatively 
flat  region,  where  the  value  of  the  dependent  variable  does  not  differ 
greatly  from  zero,  which  is  also  the  value  of  the  irregular  contour. 
Slight  smoothing  before  contouring  can  eliminate  most  of  these  irregu¬ 
larities,  as  is  shown  in  Fig.  7.  The  original  field  portrayed  in  Fig. 

6  was  smoothed  by  substituting  (8G  +  G  . ,  +G  ,  +G  + 

m,n  m+l,n  m-l,n  m,n+l 

n_p/12  for  G^  ^  at  each  interior  point  (m,n) .  The  appearance  of 

the  contours  improves  somewhat  in  the  upper  right-hand  part,  but 

smoothing  adversely  affects  the  important  small-scale  patterns  in  the 

_2 

lower  left.  For  example,  the  minimum  of  -3.87  X  10  on  the  left-hand 

_2 

boundary  of  Fig.  6  becomes  a  minimum  of  -2.81  X  10  in  Fig.  7,  and 


.  6  —  Example  of  a  slightly  "rough"  field 
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the  contour  of  value  -3.125  *  10  surrounding  it  vanishes.  Hence  this 
type  of  smoothing  before  contouring  is  not  recommended  for  any  field 
with  small-scale  features  of  Interest. 

Figure  8  depicts  a  field  with  no  apparent  irregularities.  Over 
most  of  the  area  the  value  of  the  dependent  variable  is  exactly  zero, 
which  is  also  the  value  of  the  outermost  contour  shown. 

The  field  depicted  by  Fig.  9  is  a  little  more  complicated,  but 
still  smooth.  There  are  three  zero  contours  separating  the  two  posi¬ 
tive  and  two  negative  regions,  Identified  by  the  symbols  for  maxima 
and  minima. 

The  use  of  straight-line  segments  in  these  examples  occasionally 
produces  undesirable  results,  but  on  the  whole  the  charts  are  acceptable. 
Certainly  charts  produced  this  way  can  be  of  great  value  in  enabling 
one  to  study  many  fields  without  the  slow  and  tedious  work  of  trans¬ 
cribing  and  analyzing  them  by  hand. 
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